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METHOD FOR DESIGNING THE SPIDER SPRING OF A CONTROL CLUSTER 
OF A NUCLEAR FUEL ASSEMBLY, A CORRESPONDING SYSTEM, COMPUTER 

PROGRAM AND PRODUCT 



FIELD OF INVENTION 

The present invention relates to a method for designing a 
nuclear fuel assembly which is intended to be positioned in 
a nuclear reactor, the assembly comprising a plurality of 
5 guide tubes, and a control cluster which itself comprises a 
plurality of control rods which are received in the guide 
tubes and a support for control rods, the assembly 
comprising a helical spring for damping the impact of the 
support against an upper end piece of the assembly in the 
10 event of the control cluster falling during a shutdown of 
the nuclear reactor. 

BACKGROUND OF THE INVENTION 

It will be appreciated that nuclear fuel assemblies must be 
dependable in order to allow reliable operation of nuclear 
15 reactors. Thus, design and construction provisions for such 
assemblies have been drawn up. 

These provisions impose a general framework and minimum 
criteria which the assembly constructors must take into 
consideration. 

20 As far as the helical damping spring is concerned, the 

design provisions require verification by means of tests 
that the integrity of the spring has not been affected 
during the impact brought about in the event of a shutdown 
of the reactor. 

25 Although the criterion imposed by the design provisions 
allows assemblies to be designed with satisfactory 
reliability, it would be desirable to limit the safety 
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margins during design in order to reduce the mass and the 
cost of the assemblies constructed. 

SUMMARY OF THE INVENTION 

An objective of the invention is to overcome this problem by 
5 providing a method which allows reliable nuclear fuel 
assemblies to be designed, while limiting the design 
margins . 

To this end, the invention relates to a method for designing 
a nuclear fuel assembly which is intended to be positioned 

10 in a nuclear reactor, the assembly comprising a plurality of 
guide tubes and a control cluster which itself comprises a 
plurality of control rods which are received in the guide 
tubes and a support for control rods, the assembly 
comprising a helical spring for damping the impact of the 

15 support against an upper end piece of the assembly in the 

.event of the control cluster falling during a shutdown of <. 
the nuclear reactor, wherein the method comprises the 
following steps : 

a) establishing a progression of a speed of the control 
20 cluster after the impact of the support against the upper 

end piece, 

b) establishing, based on the speed established in step a) , 
a maximum longitudinal load for compression of the spring, 
and 

25 c) establishing, based on the maximum longitudinal load for 
compression, at least a maximum shearing stress in the 
spring . 

According to specific embodiments, the method can comprise 
one or more of the following features, taken in isolation or 
30 according to all technically feasible combinations: 

- a maximum shearing stress is a shearing stress along the 
neutral axis of the spring, 
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- a maximum shearing stress is a shearing stress along the 
axis of the spring nearest the longitudinal centre axis 
thereof , 

- the method further comprises a step for verifying, using a 
5 maximum shearing stress established in step c) , that a 

maximum stress admissible by the spring has not been 
exceeded. 

The invention further relates to a system for designing a 
nuclear fuel assembly, wherein it comprises an arrangement 
10 for performing the steps of a method as defined above. 

According to a variant of the invention, the system 
comprises a computer and storage arrangement, in which at 
least a program comprising instructions for performing steps 
of the method for designing a nuclear fuel assembly is 
15 ... stored.. . :c 

The invention further relates to a computer program 
comprising instructions for performing the steps of a method 
as defined above. 

The invention also relates to a medium which can be used in 
20 a computer and on which a program as defined above is 
recorded. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be better understood from a reading of 
the description below which is given purely by way of 
25 example with reference to the appended drawings. 

Figure 1 is a schematic, perspective cutaway view of a 
nuclear fuel assembly which is designed by a method 
according to the present invention. 

Figure 2 is a schematic, partially sectioned side view drawn 
30 to an enlarged scale of the structure of the spider of the 
assembly of Figure 1 . 
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Figure 3 is a partial schematic side view of the assembly of 
Figure 1, illustrating more particularly a pair comprising a 
guide tube/control rod. 

Figure 4 is a block diagram illustrating the system for 
5 designing the assembly of Figure 1. 

Figure 5 is a flow chart illustrating successive steps of 
the design method carried out by the system of Figure 4 . 

Figure 6 is a progression curve of the falling speed of a 
control rod before it is introduced in the lower portion of 
10 the corresponding guide tube, this progression being 
calculated by the system of Figure 4 . 

Figure 7 is a progression curve of the falling speed of the 
same control rod in the lower portion of the corresponding 
guide tube, this progression being calculated by the system 
15 of Figure 4 . 

DETAILED DESCRIPTION 

Figure 1 illustrates a nuclear fuel assembly 1 which mainly 
comprises a square-based lattice 2 for nuclear fuel rods 3 
and a control cluster 4 . 

20 The assembly 1 comprises grids 5 for maintaining the rods 3, 
the grids 5 distributed over the height of the rods 3 . A 
lower end piece 6 is arranged under the lower ends of the 
rods 3 and an upper end piece 7 above the upper ends of the 
rods 3 . The upper end piece 7 is provided with springs 8 for 

25 pressing against the upper bearing plate of the reactor 
core, in which the assembly 1 is intended to be placed. 

The control cluster 4 comprises a plurality of control rods 
10, for example, 24. Conventionally, the control rods 10 
comprise a material which absorbs neutrons. 

30 The rods 3 and 10 extend in parallel with a vertical 
longitudinal direction L. 
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The rods 10 are carried at the upper ends thereof by a 
support 11 which is generally referred to as a spider . 

As illustrated more particularly in Figure 2, the spider 11 
comprises a vertical central upper head 12 and a series of 
5 arms or vanes 13 which extend radially outwards from the 

lower end of the upper head 12 as far as the radially outer 
ends 14 thereof. Each control rod 10 is connected to an arm 
13 at the upper end thereof. 

The upper head 12 of the spider 11 has a central blind hole 
10 15 which opens towards the bottom and in which a damping 
helical spring 16 is received. The spring 16 extends 
vertically along a center axis A. A tightening screw 17 
extends substantially over the entire height of the hole 15 
and is screwed into the wall 18 delimiting the upper portion 
15 of the hole 15. 

The lower portion of the screw 17 extends through the base 
. of a retaining ring 2 0 which rests on the lower end of the 
spring 16. The head 21 of the screw 17 rests, at the top, 
against the base of the retaining ring 2 0 in order to press 
20 the spring 16 against the wall 18 of the upper head 12. 

As illustrated in Figure 3 for a control rod 10, each 
control rod 10 is received in a respective guide tube 24 
which is arranged in the lattice 2 of fuel rods 3 . In this 
manner, 24 pairs comprising a guide tube/control rod are 
25 formed. Since each of these pairs has a similar structure, 
only one will be described below. 

The guide tube 24 extends from the lower end piece 6 as far 
as the upper end piece 7. The guide tube 24 comprises a 
lower portion 26 of reduced inside diameter and an upper 
30 portion 27. The lower portion 26 is connected to the lower 
end piece 6 by a collared screw 28, through which a vertical 
through-hole 2 9 extends. 
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The lower portion 2 6 of the guide tube 24 surrounds the 
control rod 10 with a radial passage gap J. 

The upper portion 27 is fixed to the upper end piece 7 and 
opens at the outside of the assembly 1. 

5 Lateral apertures 30, only one of which can be seen in 
Figure 4, are provided in the upper portion 2 7 near the 
lower portion 26. 

When the assembly 1 is placed in a nuclear reactor, the 
cooling liquid of the reactor fills the interior of the 
10 guide tube 24 . 

Conventionally, the control cluster 4 can be moved 
vertically relative to the remainder of the assembly 1 in 
order to allow adjustment of the reactivity during normal 
operation of the reactor, and therefore variations in power 
15 from zero power up to maximum output depending on the 

vertical introduction of the control rods 10 in the lattice 
2 of rods 3. The vertical displacement of the control 
cluster 24 is conventionally carried out by a drive rod 
which is connected to the upper end of the upper head 12 . 

20 When the reactor is shut down, the drive rod and the 
assembly 4 fall due to gravity. 

At the start of this falling movement, the control rods 10 
are guided only by the upper portions 2 7 of the guide tubes 
24 and have not yet reached the lower portions 26. 

25 Once the falling action has ended, the lower ends of the 

control rods 10 are introduced in the lower portions 26. The 
cooling fluid contained in the portions 26 is then violently 
forced, on the one hand, upwards thereby and, on the other 
hand, downwards through the apertures 2 9 of the collared 

30 screws 2 8 . 
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Each lower portion 2 6 therefore behaves in the manner of a 
hydraulic damper braking the falling movement of the 
corresponding control rod 10, and therefore of the assembly 
4 . 

5 This braking phase ends at the end of the travel path with 
the impact of the spider 11 against the upper end piece 7 of 
the assembly 1 . 

This impact is carried out by a retaining ring 20. During 
this impact, the spring 16 is compressed vertically in order 
10 to absorb the shock. 

According to the invention, the assembly 1 has been designed 
in order to take into consideration the specific stresses 
brought about in the assembly by the fall of the control 
cluster 4 during such a shutdown of the reactor. 

15 In. -this manner, in order to .design the assembly 1, in . .;..v 
. : particular a data-processing system 32 has been used, as . . 
illustrated schematically in Figure 4.. 

This system 32 comprises, for example, a computer or data 
processing unit 34 comprising one or more processors, a 
20 storage arrangement 36, input/output arrangement 38, and 
optionally display arrangement 40. 

Instructions which can be performed by the computer 34 are 
stored in the form of one or more programs in the storage 
arrangement 3 6 . 

25 These instructions are, for example, instructions in FORTRAN 
programming code . 

These various instructions, when they are performed by the 
computer 34, allow the method illustrated by the flow chart 
of Figure 5 to be performed. 

30 In a first step illustrated by the box 42 of this Figure, 
the computer 34 calculates, based on data 43, the 
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development of the falling speed of a control rod 10 in the 
upper portion 2 7 of the corresponding guide tube 24 in the 
event of a shutdown of the reactor. 

This calculation can be performed assuming, for example, 
5 that the control rod 10 is first subjected to constant 
loads : 

gravitational force: fg = Mg, 

Archimedes 1 thrust: fa = -p gV, 

pressure difference in the core: fc, and 

10 mechanical friction: fm, 

where M and V are the mass and the volume, respectively, of 
the assembly 4 and the drive rod thereof . 

The control rod 10 is also subjected to loads as a function,., 
of the speed or position thereof, for example, hydraulic 
15 friction which can be obtained from: fh = -ci (M + pV) v 2 , 

with v = speed of the assembly 4 and therefore of the rod 10 
in question . 

Thus, the equation of the movement of the rod in the upper 
portion 27 of the guide tube 24 is as follows: 



This gives : 

— = c2-cl v 2 
dt 

with cl = hydraulic friction in the guide tube and 
C M+pV 

25 Cl and c2 are, for example, experimental data measured 

during drop tests of the control cluster 4. These data are, 



20 
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with the other data necessary for the calculation, such as 
the mass and the volume of the assembly 4 and the drive rod 
thereof, introduced, for example, in the form of a file 43 
by way of the input/output arrangement 38. 

The computer 34 resolves the equation of the movement of the 
control rod 10, for example, using the NEWTON method. 

Thus, the progression of the speed of the control rod 10 in 
the upper portion 27 is known as a function of time. The 
profile established in this manner can be displayed in the 
form of a curve by the display arrangement 40. This curve is 
illustrated by Figure 6. 

In this manner, at the end of the step illustrated by the 
box 42, the speed of the control rod 10 is known at the 
point of entry to the lower damping portion 26 of the guide 
tube 24 . .. . 

Based on the results of the step .of , box 42 , the computer 34 . 
calculates the progression of the . speed of the control rod 
10 during its fall in the lower damping portion 26. 

This step is schematically illustrated by the box 44. 

This step can be carried out using the following equation: 



SCA = cross-section of the rod 10 and 

NCA number of rods 10 in the assembly 4 . 

Therefore, the hypothesis that f c and f m are negligible is 
applied here. 




with 



c2 = 



fg\ fa M- pV 
M + pV~ M+pV 
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The difference AP represents the elevated pressure produced 
in the cooling liquid contained in the guide tube 24, that 
is to say, the pressure thereof between the lower end of the 
rod 10 and the pressure present in the upper portion 2 7 of 
the guide tube 24 . 

AP can be established by the following formula: 
AP = — pQ 2 v 2 (EXPA + CONTRA + FECRxCISAxz) 

( SCA ( sacm\\ 2 

where EXPA-{— [l^^Jj 

with SM = cross-section of the lower portion 26, 

SACM = SM - SCA = cross -section of the annular space between 
the rod 10 and the lower portion 26, 

SACTG.= STG - SCA, where STG is the cross-section of the v 
upper .portion 27 of the, guide tube 24,. 



CONTRA = 0.4 




FECR = coefficient of loss of load owing to friction in the 
lower portion 26, 

(SCA\ 2 1 
C,SA = { ffl/J ~DM' 

DM = mean diameter of the guide tube 24 in the upper portion 
27, 

z = height of the rod 10 introduced in the lower portion 2 6 
of the guide tube 24, and 

Q = fraction of liquid, flowing upwards out of the lower 
portion 26. 

The resolution of the equations governing the movement of 
the rod 10 after entry into the lower portion 26 is carried 
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out by the computer 34, for example, using the RUNGE-KUTTA 



Thus, at the end of the step 44, the progression of the 
speed of the control rod 10 in the lower portion 2 6 of the 
5 guide tube 24 is known before the impact of the spider 11 on 
the upper end piece 7 . 

The speed profile established in this manner can be 
displayed, for example, by the arrangement 40, as 
illustrated in Figure 7. On the curve in Figure 7, the 
10 speed profile established during step 44 is the portion 
located to the left of the point 45. 

The computer 34 then performs, in the step of box 46, the 
calculation of the maximum elevated pressure produced APmax. 

This calculation can be performed, for example, based on the 
15 formula: •■ ' ' 



The computer 34 performs, in the step of box 48, the 
calculation of a circumferential stress and maximum normal 
oeMAx / to which the lower portion 2 6 of the guide tube 24 is 
20 subjected owing to the maximum elevated pressure APmax. 

This stress can be calculated based on the formula: 



where DPM = inside diameter of the lower portion 2 6 and EMP 
= minimum thickness of the wall of the lower portion 26. 

25 The system 32 can then provide, owing to the input/output 

means 38, a first result in the form of a file 49 containing 
the value Ogmax established, and optionally the maximum 
elevated pressure APmax established. 
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Next, the system 32 performs the calculation of the 
progression of the speed of the control rod 10 after it 
comes into contact with the spider 11 and the upper end 
piece 7. This calculation step is illustrated by the box 50 
5 in figure 5 . 

This calculation can be performed, for example, using the 
following equation when the ring 20, and therefore the 
spider 11, is in contact with the upper end piece 7: 

( Af + pV)^=(M- pV)g- PRCH- K(z- LAl)-c3 v 

10 with PRCH = pretension of the spring 16 = PRCMP x K, where 
PRCMP is the precompression of the spring 16 and K the 
rigidity of the spring 16, 

LAI = distance travelled by the control rod in the lower 
portion 26 before impact, and . 

15 c3 = coefficient of hydraulic damping in order to model the 
damping in the lower portion 26. 

In the event of a rebound, that is to say, when the spider 
11 is no longer in contact with the upper end piece 7, the 
equation for movement of the control rod 10 in question is 
20 written as follows: 

(M+ pV)^=(M- pV)g-c3 v 

These two equations are integrated by the computer 34, for 
example, using the RUNGE-KUTTA method. 

Therefore, the step 50 allows the kinematics of the control 
25 cluster 4 to be established during the mechanical damping of 
the shock by the spring 16. The speed profile established in 
this manner can be displayed, for example, by the 
arrangement 40. This profile corresponds to the portion 
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located to the right of the point 4 5 on the curve in Figure 
7 . 

Based on the results of this step, the system 32 performs, 
in the step 52, the calculation of a maximum vertical 
compression force Fmax, to which the spring 16 is subjected 
during the mechanical damping. 

This calculation can be performed, for example, based on the 
following formula: 



MAX 



= max (K(z-LAI) + PRCH} 



10 The system 32 then performs, in the step of box 54, the 

calculation of an approximate maximum shearing stress Tmax in 
the spring 16 : 

SF MAX DFN 
Tmax = ttDFR 3 

with DFN = DER-DFR and 
15 DER = outside diameter of the spring 16, 

DFR = diameter of the wire of the spring 16. 

Subsequently, the system 32 can optionally perform, based on 
the maximum stress Tmax/ the calculation of maximum corrected 
stresses. 

20 These stresses can be calculated by multiplying Tmax by 
different factors . 

Thus, it is possible to calculate: 

x max\ = t max X ^c' 

r MAX2 = T MAX X ^5 

0.5 

withKc = 1 + — , 

DFN 
C = — — , and 

DFR 

4C-1 0.615 
K " 4C-4 + C 
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The stress Tmaxi corresponds to the shearing stress along the 
neutral axis FN (Figure 2) of the spring 16. The stress 
corresponds to the stress along the axis F2 (Figure 2) of 
5 the spring 16 nearest the vertical center axis A of the 
spring 16 (see Figure 2) . 

At the end of this step illustrated by the box 56, the 
system 32 provides the various maximum shearing stresses 
calculated, for example, in the form of data stored in a 
10 file 57, which are transmitted by the input/output 
arrangement 38. 

Based on the data contained in the files 4 9 and 57, which 
have also been stored in the storage arrangment 36, the 
computer 34 will verify that the maximum stresses calculated 
15 are indeed acceptable for the materials' which respectively;. , 
constitute the guide tube 24 and the helical spring 16. J. 

This step has been schematically illustrated by the box '58 
in Figure 5. During such a step, the system 32 will, for 
example, verify that the maximum shearing stresses 
20 calculated during the steps 54 and 56 are less than maximum 
values admissible by the material which constitutes the 
spring 16 . This verification is performed by a comparison of 
Tmax, Tmaxi and Tmax2 with a maximum value admissible by the 
material of the spring 16. 

25 As far as the maximum circumferential stress o eM Ax is 

concerned, the verification can be performed based on a 
formula of the type: 

f (°*^MAX ) < ^admissible 

where a a dmissibie refers to the material which constitutes the 
30 lower portions 26 of the guide tubes 24. 
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The function f can be a function which takes into 
consideration other stresses to which the guide tubes 24 can 
be subjected. Such a stress can be a vertical compression 
stress o A/ to which the guide tubes 24 are subjected during 
5 the contact of the springs 8 of the upper end piece 7 
against the upper bearing plate of the core in order to 
counterbalance the hydrostatic thrust during operation. 

Thus, the function f can be, for example, in the form of f 

OA) = O0MAX + oA. 

10 It will be appreciated that this last step, illustrated by 
the box 58, can be performed by separate software which 
generally performs the validation of various design 
parameters of the assembly 1 based on results provided by 
various pieces of software each dedicated to taking into 

15 consideration specific operating conditions and which 

include the software which performs the steps 42, 44, 46, 
48, 50, 52, 54 and 56. 

In general terms, the file 43 comprising the data 43 used by 
the method for the various calculations can comprise the 
20 data of Table 1 below. 
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outside diameter of control irod 10 


(m) 


Nominal ; maximum 


inside diameter of upper portion 27 


(m) 


Nominal ; maximum 


inside diameter of lower portion 2 6 


(m) 


Nominal; maximum 


total length of lower portion 2 6 


(m) 




damping travel before impact 


(m) 




minimum thickness of wall of lower portion 26 


(m) 




maximum roughness of rod 10/ tube 24 


(m) 




diameter of aperture 2 9 


(m) 




length of aperture 2 9 


(m) 




roughness of aperture 2 9 


(m) 




mo v i ng ma s s M 


(kg) 




volume t.i~ic mass of liquid 


(kg/m 3 ) 




kinematic viscosity of liquid 


(m Ib) 




cl 


(/m) 




c2 


(m/s" 5 ) 




Young's modulus of guide tube 24 


(Pa) 




Poisson's ratio of guide tube 24 






spring precompression 16 


(m) 




preloading of spring 16 


(N) 




length of spring 16 with contiguous turns 


(m) 




outside diameter of spring 16 


(m) 




diameter of wire of spring 16 


(m) 




compression when upper head 12 is in contact 
with upper end piece 7 . 


(m) 




END 







Table 1 
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Similarly, the file 49 comprising the results from step 48 
an comprise the data of Table 2 below. 



APmax: maximum elevated pressure in lower portion 26 


(Pa) 


Zmax= corresponding penetration in lower portion 26 


(m) 


a 6M Ax: maximum stress in lower portion 26 


(Pa) 


fmax: maximum force on lower end piece 6 


(N) 


tdur: duration of fall in lower portion 2 6 before impact 


(s) 


vfin: speed of impact of assembly 4 on upper end piece 7 


(m/s) 



Table 2 

The file 57 comprising the results of step 56 can itself 
5 contain the data of Table 3 below. 



Fmax: maximum compression force on spring 16 


(N) 


• h^:' maximum deflection of spring 16 


(m) .... 

M 


t max: approximate maximum stress in spring 


(Pa) 


t*MAxi: approximate maximum stress corrected by Kc 


(Pa) . 


' t'wAxi'' approximate maximum stress corrected by K (Wahl coefficient) 


(Pa) 



Table 3 



It has been possible to verify by experiment that the 
maximum elevated pressures and the maximum stresses obtained 
by means of steps 42, 44, 46 and 48 were reliable. In this 

10 manner, the first corresponding part of the method allows 
reliable guide tubes 24 to be designed. Furthermore, this 
first part calculates only a single stress which appears to 
be the pertinent stress for the conditions being considered. 
Consequently, this first part of the method allows the 

15 security margins to be limited during design, and therefore 
assemblies which are relatively light and economical to be 
designed. 
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The second part of the method, which corresponds to steps 
50,52, 54 and 56, also allows maximum stresses to be 
reliably calculated, as confirmed by experiment. 

Thus, the second part of the method allows a reliable design 
5 to be arrived at by calculation for the spider springs 16, 
which design is found to be advantageous in comparison with 
the method of tests alone which is currently imposed by 
provisions. It will be appreciated that the second part of 
the method calculates only the small number of stresses, and 
10 in particular those on the axis F2 of the spring 16 nearest 
the center axis A of the spring, which are found to be 
pertinent to the conditions envisaged. In this manner, the 
second part of the method allows the design margins to be 
reduced. 

15 In more general terms, the steps 42, 44, 46 and 48, on the.,, 
one hand, and 50, 52, 54 and 56, on the other, can be 
performed by separate pieces of software. 

In order to. increase the reliability of the calculation, for • 
performing the first part of the method it is possible to 
20 use, as the passage gap J, the nominal value of the gap, or 
this nominal value corrected by the manufacturing tolerance, 
or a value resulting from statistical studies of the 
distribution of passage gaps J obtained in constructed 
assemblies . 

25 In a variant, it is possible to use a gap value J which is 
greater for steps 42 and 44 and a smaller gap value J for 
steps 46 and 48. This allows a high stress value Qomax to be 
calculated because the speed reached during the fall of the 
rod 10 in question is high and the volume available in the 

30 lower portion 26 for the liquid during damping is small. 
However, this high stress value is not unrealistic and 
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therefore does not lead to unjustified design margins, as 
illustrated by the following example. 

According to a specific variant, the upper value can be a 
maximum value for gap J which is verified with a given 
5 probability, for example, 95%, in constructed assemblies, 

and the lower value can be a minimum value obtained with the 
same probability. This variant allows an approximation of 
the situation where a single pair comprising a guide 
tube/control rod has minimum gap J, where the maximum stress 
10 ogmax would be reached, and where all the other pairs 

comprising a guide tube/control rod have the maximum passage 
gap J, which would be the most extreme case. 

In some variants, the first part of the method could also 
take into consideration forms of the lower damping portion 
15 2 6 which are different from. those described previously. In vk 
this manner, these lower damping portions could have a 
.9 plurality of successive portions of reduced diameter. 

£ optionally separated by portions of increased diameter, 

generally referred to as cavities. In some variants, the 
20 first part of the method is performed with collared screws 
28 which are not perforated by holes 29. 

In still more general terms, the first part and the second 
part of the design method described can be used 
independently of each other. In this manner, it is possible 
25 to perform the second part relating to the design of the 
spring 16 without referring to the calculation of the 
elevated pressure AP and the stress oeMAx- 
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